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4 • • SUMMARY 

A method for calculating the upwash -velocity component in the 
vicinity of thin wings at supersonic speeds is presented. The 
method is applied to obtain an explicit expression ‘for the upVash 
over wing tips of fairly general plan form and profile. As a 
special case, numerical values are presented for the slopes of the 
streamlines off the tip of a rectangular plan-form thin flat- 
plate wing. The formulation is extended to give a method for 
obtaining the velocity potential at points on arbitrary thin wings 
influenced by isolated or interacting external flow fields off the 
wing plan form, . The solutions, obtained by the method for regions 
influenced by so-called subsonic trailing edges do not conform, 
however, to the Kutta-Joukawski condition. 

The. method in principle may be applied to obtain the- aero- 
dynamic coefficients and hence the lift distribution for thin 
wings of arbitrary plan form. and profile; the calculus involved 
in obtaining explicit solutions, however, is likely to be diffi- 
cult and impractical. The functions were therefore altered to 
Isolate and to remove nonessential singularities. The equations 
so obtained ere suitable for numerical calculations of the aero- 
dynamic coefficients of arbitrary thin wings at supersonic' speeds. 
As an example of . the method, the upwash between the leading edges 
and the foremost Mach waves in the plane of the flat -plate delta . 
wing was calculated and compared with approximate results obtained 
by neglecting flow-field interaction. The pressure coefficient 
on the surface of the wing was likewise numerically computed and 
compared with the exact solution,- 

c 

“ ' INTROnJCTION ■ • ‘ ■ 

A method for obtaining the lift, drag, and pressure distrl- . 
butions in the vicinity of thin wing tips at supersonic speeds is 
presented in reference 1. The basis of the method was -to plaoe a 
thin' diaphragn along a stream sheet in the plane of the wing 
between the wing boundary and the foremost Mach wave. In this 
manner an integral equation was established to define the slopes 
of the streamlines in the disturbed flow field. The interaction 
affects of the two surfaces of the wing in the vicinity of the 
tip were thus evaluated. 
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The main object of reference 1 was to present a method for 
calculating the theoretical aerodynamic coefficients on the sur- 
face of a wing. These coefficients were obtained .without an 
explicit solution for the slope cf the stream sheet for points on 
the wing surface influenced hy noninteracting external flow fields. 
If the slopes of the streamlines could be determined, the upwash- 
velooity components in the external flow fields and the aerodjrnamic 
effects of interacting external flow fields could be evaluated. 
Pressure distributions and the lift and drag coefficients could 
then be calculated for arbitrary thin wings at supersonic speeds. 

The present report shows that the defining equation for the 
slopes of the streamlines off the wing -plan -form boundaries is a 
special case of Abel’s integral equation. The solution to Abel’s 
equation is applied in order to obtain the slopes - of tb© stream- 
lines off wings of fairly general- plan form, and profile. As a 
special case, numerical values are presented for the slopes of the 
streamlines off the wing tip of a rectangular plan-form thin flat 
plate. The general formulation is shown to give a method for 
obtaining the veloaity potential of finite wings influenced by 
interacting external flow fields. The solutions obtained by the 
method for regions influenced by so-called subsonic trailing- edges 
do not conform, however, to .the -Kutta-Joukowski condition, 

Although the method in principle may be applied to obtain the 
aerodynamic coefficients and hence the lift distribution for 'thin 
wings of arbitrary plan form and profiles, the calculus involved 
in obtaining explicit solutions is likely to be difficult and 
impractical. The functions are therefore altered to isolate and 
remove nonessential singularities. The equations so obtained are 
suitable for numerical calculations of the aerodynamic coeffi- 
cients of arbitrary thin wings at supersonic speeds. As an 
example of the method, the pressure coefficient of the delta wing 
included within the Mach cone is computed by Numerical methods and 
compared with the solution "obtained by other methods . 


SLOPES OE THE STREAMLINES 

The analysis is considerably simplified in a set of oblique 
coordinates u,v whose axes lie parallel to the Mach lines of the 
flow, 'In this set, the value -of one of the coordinates of a point 
is the distance measured parallel to the coordinate axis from the 
point to the other coordinate axis. If the wing lies in the 
x,y plane and the free-stream velocity is parallel to the x axis, 
the transformation equations from one set of coordinates to the 
other are 
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u = |r ^ 7 v= ^ (* + py) 

* = |' (u - + v) y = H (v “ u) (D 

where M is the Mach number and- 0 = J$-l. (A list of symbols 
is included in appendix A. ) 

If the coordinate zero of the two systems, is plated on the 
point of tangency of the ying boundary (fig. i) and the foremost 
Mach. wave and if the wing boundary is represented by the two 
equations v = v 1 (u) and v ?= VgCu), the defining equation for. 

the slope of the stream sheet (referred to as a diaphragm) 
near the plane of the wing measured in y = constant planes is 
given in reference 1 as 



where u D and v D represent idle coordinates of a local point on 

the diaphragm and op and Orj are the slopes (measured in 

radians) of the wing on the bottom and top surfaces, respectively. 
(The sigh of the wing slope is defined oppositely on the bottom 
and top surfaces; for example, c- g and Orp are both positive for 

a wedge-profile wing at zero angle of attack.) 

The derivation of equation (2) presented in reference 1 
required that the three components of the perturbation velocity 
be continuous across the diaphragm. . If the wing sheds a vortex 
sheet, discontinuities in one or more of the perturbation-velocity 
components are feasible. Because the diaphragm can sustain no 
pressure difference between its. top and' bottom surfaces, the 
x component of the perturbation velocity must be continuous. The 
partial derivative with respect to x of equation (2) then applies 
rather than equation (2). Integration yields equation (2), except 
that an arbitrary function of y may be added to either member. 
This funotiop of y represents the circulation in. the yortex 
sheet and may be adjusted to make the solution to the equation 
satisfy the Kutta-Joukowski condition in cases involving so-called 
subsonic trailing edges. 
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Inasmuch, as the limits of integration with respect to u are 
the same for "both members of equation (2) for all values of Up 

and because of the nature of the functions, the equation may he 
reduced, as in reference 1, to the form 



(Cj-OrjOdv - 


rj(u) ' 2 *Ad-^ 


(3) 


Equation (3) is a special case of Abel's integral equation. 
In the; notation of reference. 2, if 


■'X 0 < fi <; 1 

d 4 = f{x). . . 

a (x-g) 1 * a x b 


(4) 


where Xja, and b. are real and finite, then "the continuous 
solution, if It exists, can be- none other than" 


u(z) a !^4l 

it dz 


f (x)dy 
1-u 

(z-x) 


. The following table compares the notation of the symbols of 
equations (3) and (4): ' 


Notation of 
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(The function f ( Vp) -is usually not zero at .v-p = v 2 {u) nor is 
f 1 (irjj) ■ continuous^ as - required in reference '2. A singularity 

exists, in the present problem that maybe isolated. For a discus*-' 
sion of the requirements, for application -of -equation {5) to obtain 
a solution of equation (4), see reference 3.) ‘The- solution “af 
equation (3) isi then given by equation (5) as ‘ • •• ■ ' 


* (u * z) " §3? 


dv- 


D- 


,v 2 (u) 




rjzbr, .. 

' ‘ ( < ^B~0»ji)dv 

, v 2 >d ~ T- 
T l (u) • :: ■■ , 


( 6 ) 


evaluated at z = v. 

The upwash 
by the equation 


The upwash -velocity component V_ id 'directly related to A 

4 • i 


v 2 * XU. 


(7)' 


where U is the free-stream velocity.. Inasmuch as the sign of A 
vds chosen Tflth respect to the top surface , of the diaphragm, a 
positive A implies an upflow. Substitution of A from equa- 
tion (6) into equation (7) yields ! 



~z 

rr 2 (u) . 



dT D 


(8) 

z n 5z 

w 

v,(u) ^ J 

2yv B -T 

^(u) 


evaluated at, z » v. 


Equation (6) is relatively easy to evaluate when the profile 
of the wing is symmetrical about the plane of the wing. In this 
case, (ap-ap) - 2a' where a is the angle- of attack. The 

indicated integration for the wing plan form of figure 1 is car- 
ried out in' appendix B as equation (B5j to give 


A(u,v) 



/V 2 (u)-V;l(u) 
'' .V-Vg(u) ' 


r tan 


■1 /vg (u) -v 1 (u) 

Y J ,T-Vg(u)" 
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That equation (9) is the solution of equation (3) for the symmet- 
rical wing has heen verified hy direct substitution. (See 
appendix B.) -(If o B - c T is a function only of u, the solution, 

. - * t v 

equation 9, for ■ A still, applies providing the factor. 2 a is 
replaced by 03 - a^. ) • 

A plot of A /a (from equation (9)) as a function of the 
Yofuj-v-i (u) 

independent variable — is presented in figure 2. Along 

« . ' VrV 2 (U'J 

the foremost Maoh wave] v 2 (u). - v 1 (u) = 0 and X = 0 . The value 

of X becomes infinite at the wing boundary corresponding to the 
mathematical discontinuity along V = v 2 (u). 

The quantity may be interpreted geometrically 

vtv 2 (u) 

with the aid of figure 1. The quantity v 2 (u)-v-j_(u) is the length 
of the line a b, whereas v-v 2 (u) is the length of the line b v. 
From figure 2 it is then apparent -that v must be close to the 
wing boundary before A/a becomes very large. 


may be interpreted geometrically 


As a further .illustration, the diaphragm slope of a rectan- 
gular plan-form thin flat-plate wing has been calculated. Id this 
case, v-j/u) = -u, v 2 (u) = u, and by equations ( 1 ) 

^2^) ■- 71 ( u ) 2u x , 


v-vg(u; v-u py v ' 

The quantity py/x was therefore 'taken as the independent variable. 
The resulting slope ratio A/a is presented in figure 3 as a 
funotion of py/x, The ratio A/a remains fairly small except 
near the wing edge. For example, at a Mach number of 3, A /a is 
still only 4.26 for y/x = 0.0053. 


APPLICATIONS TO AERODYNAMIC THEORY 

The explioit solution (equation ( 6 )) of the integral equa- 
tion (3) has direct application in the calculation of the lift 
distributions of thin wings. The solutions obtained for cases 
involving so-called subsonic trailing edges may violate the Kutta- 
Joukowskl condition and hence should be used with caution. In 
reference 1 a method -is presented for obtaining the' ■aerodynamic 
coefficients in the vicinity of wing tips under the influence of 
independently disturbed external flow fields; that is, each dis- 
turbed external flow field includes no other external unknown flow 
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field la -lts.forwapd'-Mach cone. By application' of equation (6), 
portions of the external.- flow field may he* evaluated, The- inter- 
action effects of one external field on another may he determined# 

The details of- the 'process' frill he illustrated' for the. wing 
plan fora of figure d. The. leading edgea are defined : by, the_- ■ • 
equations’- 


v • 

u 

> 

or u = u^(v) 

V a V 2 (u) •• 

or u ® ug(v) 

V s v 3 (u) 

or u * 113(7) 




The diaphragm area Sp may he divided into sections Sp and . 
®D,2 j where there’ Is no interaction,, - and Sjj 3 and Sp where 
there is interaction.' The.. slope N of the diaphragm Sp g is 
given hy equation (6) as . - 


f*z 


>!><»,*) 


dvr 


Jy 2 (»>. 


a/2-VD 


• • - (Pb^T ^ 7 


evaluated at z ' «= t T In a similar manner, the slope of the 
diaphragm S-^ j is ; . • - " 


- 1 g- 


du D 


> 3 (V) 


u 3 (v) 


U X (v) 


(? B -q r )du 

2 j/u-j) r h 


( 11 ) 


5, ir 


(12) 


evaluated at z » u. 


Either of two schemes may now he applied to extend the calcu- 
lation of the velocity potential beyond the shaded region. 1 / In 
calculating the influence of., the external field ‘"Sp g and Sp 4* 

the diaphragm S D ^ is .considered as part of the wing; similarly, 
when the' influence 1 of the external 'fields - and Sj^-g is tp ? 

be calculated, the diaphragm £§p |2 1® considered -as part, of the-; 
wing. The methods of reference l,.then directly apply. 

** v> - • ’ 

•w • 

* * « * 


■» 



, velocity potehtlal-^tjf at stay point- * v w . : . ^the top 

sur^eio& at -the vine is given 'by the' integral - - «,. 

..S': . . * .1 , ■ 


•- 1 




t ‘ * 

u , .;Cf^du ; dy 


rr» 


*w 


y (Uw-u) (V w -v) 


u_ 

jtM 


Xdu dv 


% 


V(vtt)(vv) 


(is) 


or in Cartesian coordinates 

i " ■ ■ * I, 




*'£ •- ■'~ : U 1 


• ^ oi?c, g^- dl’d-'n - n -• 

w fry. . • ‘ < . 

®w * % 


• % n-.'Xdfcidb'H 

|%>. •, 

. * *• 1 - *-*>'{.•• • . . 

(14) 

■where TJ is the free-stream velocity (parallel to the x axis) 
and and- Sjj are the wing and diaphragn. areas included in the 
'forward Mach cone from the po'ipt v w ', or. x,y.- 

‘ r'S ■* ' ; 

If for simplicity the quantity E? .’is defined as 

• ! *iS .*. • v; :v. v>s if.-.lv %r nl - !■ :: c . . . . . 

*(x,jr,|,i|) s 


n/ (x-|T 2 ^ 2 (y-n) 2 £ ‘‘ " 


1 •» ^ 

then in the.^t.ati.^i hr figure 5, 

* ■- i . r . 


i i. 1 


rr> 


q> T = - 




S^l+2+3+4) 


a^dt drj - 




Spd+3+3+,4) 


XKd£ dq (15) 


■'* :j ^ V r v _ . 

where the notation ^^(i+2+3+4) --means integratiOp. over the sur- 
face- S^ i-, + Sjr 2 .+; t « . .According .to 'the. methods, of refer* 
©no© 1' - - . 


s 

« • ► « -L 

t . 

*s \ ‘ i 

[>>* 


“ Tv'E-d^'dii =:■ 

X ' 

Jw 

* * * m * •) 1* » ‘'AJlf* 1 i. v»« • f A 

• .■ v »•* \ ■> : ; • •• 

Sj( 1+2+4) 

T .* 

' £ 

V 


in (16) 


l^(US) 
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(Cfp-Om) 

. - - .AEd| cLtj S - -- -- - Kd£.dq 

J^S^ 1+2+3) jjS^Cl+2) . 

Addition of equations (16) and (17) yields 


AKdtdi} ~ 

l&p(l+2+3+4) ^ S^r (2+3) 


. (CB'-Orp) _ t , 
Kd^dtj- + 


Sw,l 


; . AKdldq 
£^(1*2) • 


Substitution of equation (18) into equation .(15) yields 


qp T = - OgKdfdTi - 

v v 1 ^ «> ®w, (2+3) 


(OB+Oij,) 


Ed|dq 


a jEdldr) + AjKd^dt] + AgEdldq 

^,4 JJ S D,1 JJ S D,2 


where A;l and Ag aye given by equations (12) and (11), 
respectively. 

An alternate scheme, would be to calculate the slope pf the 
diaphragm in regions S^^ and -S^ 4 t By applying equation (6), 
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^•i'k 


v 2 (u) 

fV u) 


fr*p- 

a/.Z" v P 


P T 3^ ^(UjVjdv . 
• ■ 


(c B - 0 T )dv 


J v 3( u ) 

evaluated at z-..= V, Similarly 




( 20 ) 



•'Z 

A 




i/z-UD. 


* 

Ugfv) 

Jo 


ugM 


-?^(u,v)du 

-/u D -u 


0 / . 


pu 3 (v) . 


+ 
r~sj 


u 2 (v)- 


(a B ~a T )du 

2 - u 


( 21 ) 


evaluated at z = u.. This process may "be continued until the 
slopes of all the diaphragms are Jmown.. The velocity potential cprp 

at any point (u w ,v w ) on the top surface of the wing is then given 

"by equations (13) and (14). 


One f urther case will he mentioned. If v = v 1 (u) reduces 

to a point at the origin, there is a continuous interaction of the 
two eternal fields. (See fig. 6.) The slope of the diaphragm in 
the region Sjj 2 may he written as 


r z 


A 2 (u,z)‘= i A 


<lv. 


D 






rV u) 


-A^(u,v')dv 


Jo 


rv 2 o*) 


Jv 3 (u) 


(c^-q T )dv 
2 t/vjpr 


(22a) 


evaluated at z ?= v or 
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A 2 (u,z) . ii 


, JV 2 {u) 


Z : 


dv T 


1 iL 

* dz 




dv T 




rv*> 

(g B -gj) dT 
.(u) 2 Sy* 

O V . . 


evaluated at 2 = y. Similarly, ' ' '■ '• * 

‘u 3 (v) 




, &*p. . 

.bm v 5 ^ 


(Og-OjjOdu 

vi) 2 ^i. 


(22b) ■■ 


1 JL 

* dz 


J5 D_ 

^(v) 


A g (u,v)du 


-/V"* 


(23) 


evaluated at z =» u, 


Thus there are two equations for the two unknown, functions A-^ 
and Ag. The functions A^ and Ag can therefore be determined, 

at least in principle ^ either by successive approximations or by 
direct substitution and solution of the resulting integral equa- 
tion. (XCTthe wing has a symmetrical plan fi^rm about the x axis, 
h 2 (u,v) = A^('v,u) and the two equations ‘(22) and (23) become a 

single integral equation. Furthermore, if the flow is conical, 
then Ag is a function of v/u. ) 


NUMERICAL CALCULATION OF DIAPHRAGM SLOPES 
AND OF PERTURBATION VELOCITIES 

. Calculation of diaphragm slopes , - The calculation of the^dla- 
phragk slope A and of the^'pertinhation velocities R? and 22 

3x 9y 
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in closed form or in series will generally ."be difficult for a wing 
of arbitrary form. Beoause the -integrals’' defining X and cp 
have finite Regions "of- integration they are well suited to purely 
numerical integration*. * * : ■ : • .. . 

In order to carry out numerical integration, it is always 
neoessary to eliminate infinite singularities from the integrands. 
For the calculation of perturbation velocities, it is furthermore 
desirable to eliminate the "extra step of the, numerical differentia* 
tion of an integral, r* ' " 

The evaluation of the integrals Involving Singularities is 
disoussed in appendix B, Equations are derived for the calcula- 
tion of X and of the derivatives of cp that are suitable for 
calculation with adding calculators or with "the aid of mechanical'- • 
computors. 

” j * 

The subsequent "equations have been derived for a general wing 
with a mean camber Iljn© that is curved at the edges of the wing. 

If the curvature of 'the camber line is zfercf at the edges, the 
terms involving the derivatives of (a B *’0’ T ) become zero. In the 

case of an uncambered wing, • (ob-^t) has the constant value of 2a. 
In these important cases the. equations are simplified. 

The final equations for :Xg(u,v;) equivalent to equations (20) 

to (22a) are derived in appendix C as equations (C13), (016), (024), 
and (C25). (See figs. 5 and 6.) 


Nj(u,v)’i - 


1 .32.0^2) _ l- r . „ , ... 

* ■"/—F 1 : 2 ••• 




■ r 


P.v • 

h* (vjjJvh? (v)’ 


it : 

.* o 




^(v) (24) . 


where 









+ A, 

*1 t d -v 3 


( A 1 V a-v 
^ 1 v^y-t^) 32 


- y^ 


fl-Cag-Oj) 

dr 


(gg’-gj) •*{<^B‘<>p)g~(Tp-72) 

4(Td- t ^ 2 . 


dv J2 


y^ + jvZ 

g2(u,T 2 ) = (cJ B TO T ) 2 ; .A 2 -V3 - Aj_ log^ 


V JJ V O 


-ft- 


Ai ^ dv + . j' * ■ 


» u. » • 

v.fe V u > ; 




and frr^-fTmln ma wriH fhft Trains of fir «+. mrmrs'p /Roo ■4Mr» A 1 

*'jj ~x'c — '"JJ ~"X' — t _ — * \ ---o* ■"•/ 

Calculation of surface velocities. - The veiooity potential cp is given In u-v coordinates 
"by- the double integral 


i 
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cp. -JL 

^ • jm 


>u w 


v w 


adu dv 


0 ,/(iv-u)(y w -v)' 


(28) 


where u^v.^ are the coordinates of a point on the wing and where 

G is to he interpreted as the wing slope, or as the diaphragm 
slope, as required. 

The velocities Y x and Yy. ' are given hy 


'.-g 


i 

_ 9qp M (b<V Sep \ 

= Sy = 2* A^w " S’W 


(29) 


The partial derivative .• may he confuted hy use of. the folloWr ; 

° u w • 

ing equation^, which are derived in appendix C as equations (033) 
to (036), . : 


dep _ U 

. §u~ 2Mjt 


lU w 


. 2R(u„,v w ) 
du - 


0 v -») 3/8 




(SQ) 


where, for points on the top surface of the wing, 

+ Jvz 

BtK^w) = 2a T (u,v w ) ,/v w -v 3 + A x (u) log e 


rv. 


go 

pv. 


k . _ii_\ 

\ 1 y v 3 -v/ -Aw^ 


v o T (u,v)-<J T (u,v w ) 


v 3 


■FF 


dv 


(31) 


r* 
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for points on the "bottom surface of the wing, . 


V^V = 2^(u,v w )V^rt ---^C u ) lo S, 




nv. 


rv. 


X _ ^1 \ , dv - + 

,1 ~ '-/V 7 ' 


w 


?3 ' ■ 


ctb(u, T )^(u; t ^ av (32) 

. r*r . 1 

/V i-T • ■ 

V .V • 


For points, on the after (.diaphragm (Sj, 4. to fig* 4),' .. Y ’ ' • 

%(-u,y D ) » - %(u,y d ) « (a B -cr T ) 2 ^-Tg - 

“ v 2 . . . „ . 

- + ^-0*) - (%-a ? ) 2 ] •' (33) 

■ . w 3 

and, for points on the forward diaphragn (S^ j In f ig.- 4), - B = 0. 

The function B is a first integral' of the equation' for the 
perturbation. potential or for its derivative -with respect to .u^. . 

A similar function may be defined by interchanging the roles' of - 
u and v. in equations (31) to (33)' to obtain .either the potential . 
or its derivative with respect to y by a subsequent integration, 

•* • ’ ‘ „ . Jv . * **■'•' 

^cp 

The paths of integration used in the calculation, cf 2~-t .are . 

• “ d“w 

shown in figure. 7 fr Path 1 is used to oompute the function B; 
path’ 2 , to compute ’ cp or 2J_ , ' . . 

... ' * 

It is noted that B 'is singular at the'point v = v 3j '' if ■. I’ . 

Aq ^ 0. This point does not appear in the integral for cp. or ■ 
^cp 

— . except when y ?= v 3 (u) and the. singularity therefore intro- 

• • ■*. . - '■* 59 

duces no difficulty, because at v = v*(u) cp - is finite and' • • 

- , - . ........ •• • - .clV- . 

is infinite to A^(u) ^0.-. ■■ - ... .• .. .■ 


j 



The integrand of equation (30) is infinite at u p u^. . It is 
not convenient to remove this singularity^ became. .to do. so would 

require the- calculation of Details of the results of 

. . du v • ; 

integration by the use of a power-series expansion of 

R(u,y w ) - R(u^,v v ) in a small region of . u . near u^ 'are given 

in appendix C -in a form suitable' for numerical "calculation. 

As an illustration of the results of the numerical method, 
the diaphragm slopes and the pressure coefficients were computed 
for an unyawed flat-plate delta wing included. within the Mach cone 
from the vertex: The locus of forward edges of the wings are 
specified in the oblique coordinates u and by the relations 

Vg = ku 

■ .Vg = u/k or Ug *t kv 

In this case the diaphragm slopes are always subject to inter- 
action effects . The calculations were carried Out by a method of 
successive approximation , Eight equally spaced stations were 
seleoted from the forward Mach line to the edge .of the wing along 
a l ine u = constant on the diaphragm. Integrations were car- 
ried out by Simpson's rule. 

t The diaphraga elopes for k » 0.25, 0,50, and 0,75 are 

shown in figure- a together with the slopes neglecting all inter- .. 

action as computed by means of equatio n (9). The dimensionless 

slppes . 'hf ct have been multiplied by* A - Xl in order to avoid 

.. v ini 

representation of the pole at ■ _L = -1. The intercept of the curve 

ku ^ 

at v/ku .= 1 is .the value of ,. -A. ; - ; /. 

The results of this integration show that: even, for a sharply • 
pointed wing (k » 0.75) the effect of the interaction is quit® 
Small. For many calculations it will therefore be permissible to 
neglect, at least for a first approximation, the -effect of the 
interaction of external flow fields' -in the. calculation of dia- • 
phraga slopes and pressure coefficients . • 

The pressure coefficients, were., canqjuted'. by the use qf values ■ 
of R ' at eight' equally spaced points’ 'along' a line v = constant, 
with two additional points plaqed near the edge... v «?. u/k. -The 
exact pressure coefficients as -computed' by the methods of refer- 
ences 4 and 5 are compared with the numerically computed pressure 
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coefficients in figure 9. Agreement', "between the pressure coeffi- 
cients as computed "by the two methods is close, even though only a 
few points were used in the calculation. 


Flight Propulsion Research Laboratory,. 

Rational Advisory Committee for Aeronautics, 
Cleveland, Ohio, August 25, 1947. 
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.SYMBOLS 


The following symbols are used in this report: 
A, B coefficients • 

a. b point values of function 


g(u,v) 

g« (u,v) 3 


pressure coefficient 

first integral for calculated diaphragm slopes 


h’ (u,v) 


regular terms of g* (u,v) 


constant > 0 


B(u,v) 


Mach number 

first integral for calculation of perturbation 
potential or its derivatives 


u, v 


x, y, z 


plan-form area 
free -stream velocity- 

oblique coordinates whose a^es lie parallel to 
Mach lines 

perturbation velocity 

Cartesian coordinates (also used as subscripts to 
indicate components of velocity along coordinate 
axes ) 

angle of attack 


p 3 /m 2 -1 cotangent of free-stream Mach angle 
8 interval of variable of integration 


Cartesian coordinates 
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X 

slope af stream sliest near plane of wing measured in 
y = constant planes 

A 

strength of singularity in diaphragm slope at •wing 
edge - • 

a 

slope of wing surface with respect to x,y plane as 
measured in y = constant planes 

qp 

perturbation velocity potential 

f 

Subscripts : 


B 

bottom wing surface 

D 

diaphragm 

T 

top wing surface 

W 

wing 

1^ 2^ 3 ; ©to. 

refers to numbered curves or surface areas 

Examples j 

. 


slope of top surface of wing 


potential on top surface of wing due to diaphragm 

^W,5 

wing area 3 . 

V 1 

curve v = v-^u) 

*1 

curve u *= UjXv) 


slope of' diaphragm in plan area 1 

(°j3'* a T)2 

i 

difference between slopes of 1 bottom and top wing 

surface at curve 2 * ■ 

» 

S2,d( v d) 

portion of function g 2 (v D ) due to diaphragm 

elements lying ahead of the wing on a line 
u p constant 

%Kv ¥ ) 

first integral of perturbation potential for 
points -on wing having the coordinate v w 
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appendix b- 


’minATION^QF ELATION (6) FOR -WING WHOSE PROFILE 


IS SYMMEERICAL ABOUT PLANE OF WING 


In this case, 
tion (6) ‘becomes 


°B~°T1 


= a, the angle of attack'; thus eq.ua- 


A(u,z) *= J 

, 11 bz 


v„(u) 


v 2 (u) D Jv-^u) 


a dv 


Evaluation of the inside integral yields 


A(u,z) 


2a d 


'v-p-Tg - 


Jv 2 (u) 


The two integrations indicated in equation (B2) may he obtained * 
from formulas 111 and 113 of reference 8 to give 


- P z 


'V V 2 


dv D = £ (z-v 2 ) 


T P (u) 


—E==k dv = yCz-ViHVg-V!) + ’(z-v-l) tan" 1 (B4 

yz-TD ... V v 2. 1 


k 2 (u) 


Substitution of equations (B3) and (B4) into equation (B2) yields 


A(u,z) 


2a 

it 


a L 

57 Is 


I (s-V 


y(z-v 2 ) (T2-T!) - (z-Vi) tan -1 ^ 


p^l 

\/ v 0 - V ^ ) 
V & XJ 


from -sfliloh 


_ / 




r= — =— \ ■ 


^ ■ - (~ i 4% '<(%-*)■ v^ij 

That equation (B5) 5s a solution of equation. (3) may "be verified by substitution: 


(B5) 


f* V 

*D 

P V D . ‘ 

r*p 

■%Vd 

dv 

X(u,v)dv ?qt 

• .1 ^-1 /. T_T 2 1t ?a 

i ^ ^ 

1 1 4 = 

/?—* — Jt 

. . -AtT** • 

. . Arr* ■ * l 

v(v-v 2 )(vd-v) 

•An-v 

t 2 (U) J 

*2^ : J 

*2< u ) 

^2( u ) 


(B6) 


The first integral of the seoond member may he integrated “by parts to give 


V D . 


1 . -1 I< t - v 2> 

■■■■- tan a ~ — ~ — 
/vjpv / V T 2~ v l 


dv » -2 yVjj-v tan' 


.1 j v-v 2 ‘ 


• it 71 *. 

riJ 


-1*2 


( V - Tl ) 


dv 


^ /(vp-rKvg^) 
tv^) yv-v 2 


dv 


(B7) 


to 

H 


OH NJj VC5VH 
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Combination of the second integral of equations (B6) and (B7) 
yields 


2a 

ir 


'D 


Jv 9 (u) 


i , -i r ™2 r 

V v 2-^i 




v 2 )(tD“ v ) 


dv 


= T v^-i) <Ti" T i } 


- _ • * 

r- 


dv 


(B8) 


(V-T a ) ,/ -y*+ (y D +y 2 )v-v 2 v D 
Equation *(B8) was integrated "by formula 195 of reference 6 to give 


P V D c- 


2a 

it 


v 2 (u) 


— i— tan -1 /-^+ /t T . dv = 

La/ J (B9) 


Evaluation of the third integral of equation (B6) yields 


pv T 


- a 


Ivo(u) 


- *» ~ 2a i/vp-Tg 


(BIO) 


Substitution of equations (B9) and (BIO) into equation (B6) yields 


M^v)..fl c 2a 

t \ . / ( v-r^-v\ N ' 


(Bll) 


>(u) 


But 


'Vp(u) 


ta (u) 


qdv 


= 2a 


j^V T D~ v l '■•:V v D" v 2^ 


(B12) 
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A comparison of equations (Bll) .anct '(’B12) shows that A as given 
"by equation (B5) ^ is a solution' of the equation 



J 
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APPENDIX C. 


NUMERICAL CALCULATION OF DIAPHRAGM SLOPES AND VELOCITIES 

In order to compute diaphragm slopes and velocities or pres>- 
Bure coefficients by numerical methods, . all the infinite singular- 
ities -with the integrals must "be removed and evaluated separately. 
Whenever possible, numerical differentiation should be eliminated. 

It is also possible to evaluate integrals having infinite singular- 
ities by expanding the integrand in series near the pole and then 
evaluating the integral by term-vise integration of the function. 

The first method has been chosen because in many cases, for 
example, for wings with flat surfaces, the integrand left after 
removal of the pole vanishes for all or part of the region of 
integration. 

The first part of this discussion in append!* C is devoted to 
the isolation of singularities in the integrals defining A. Hie 
second part considers methods of computing velocities on the wing 
surface. 

Calculation of A. - The diaphragm slope X is in general 
given by an equation Of the form of equation (20) or equation (22a), 
which for the region aft a wing surface in the direction of increas- 
ing v is 



->z 


"V 2 (u) 


-V 3 (u) 

i jL 

dT D 


a B~°T 

dy - 

A x (u,v)dv 

* dz 

% y z " v D 


/ t D" t 


V 

v 2 (u) 


k 3 (u) 


0 

/Z 


(Cl) 


evaluated at z = v, where A-^ is the diaphragm slope ahead of 
the wing on a line u = constant. 

Equation (Cl) may be formally simplified by defining the 
function ggfap) 


pv 2 (u) 


s 2 (V = 


v 3 (u) 




dv - 


■V 3 (u) 


A^dv 


/V 


( 02 ) 
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•The . calculation -.of • Ag ■'■•'majr be divided ^into two* -parts : the 
calculation of A-g^-'- ! due J -to‘-the--f tulctioh > .• ' • ' 

, s , 

■'(a B -a T ) dv 


■;S w- 


Mu) 


* 

and the calculation of A 2 p 1 due 'to the function- 

1 . « . 

p-v 3 (u) * ' ' 1 

•- «A, ~dv 


(03) 


S 2t D( v u) * 


£04) 


r 


Each of the integrands of equations (03) and (C4) may contain 
poles of order 1/2 for some value of v. The integrand defining 
32,\t- i£J .singular for v D =. v 2t (u),. .. The effect of the singularity 
may ^e isolated and evaluated by - •' •’ ' \ 1 ‘ 


S2^ W (V = 






Jv 3 (u) 


V 


V D ^V 


-Yg(u) '• . •; . A . .. 

■ " " • fa^: dT 

, , ^ 


» » ' ■ # 5 • * 


4 * 




rM u ) 

‘ ^j’ d7 v ■ 7 * 1 ‘ ,; ' 


k(u) 


2 /5r,. 


(C5> . 


t ■ 


whepe thd'tefm : '- is ' used ; to mean the valu'd ’ -A^B^rp) at 

curve ' 2 'and'"in Tfhich ■tee' second* : 'te±w has : a finite integrand -even ■’ "■ 
at Vjj = v&. 1 - * - *••• - ■ ■ ■ =• 


■i . •• 
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The integrand defining - gg jj(tjj) ha& a pole due to A^ of 
order l/2' at v D = Vg(u). The effect of this Singularity niay be 
isolkted by 


V T) 


Ai • 

yVgM-V 




u(v)- • 


(C6) 


and "by rewriting equation (C4) in the form. 


pv 3 (u) 


S2,X)( v d) 


-A^ dy 


pv,(u) 


JO 


7[js(")- T ](y T ) Jo -/V v 


-fefofe (C7) 


If the value of la computed by the use of the equation 


4l " 


lim 

v-»v 3 (u) 


> 3 (u)-V ^ 


v < v-^u) (C8) 


the funotion, |i(v) will be everywhere finite, as may be seen for a 
special case from equation (9). The function A-^ is defined by 
the equation 


rsz 


‘1 * 


1 a 


dz 




^(V) 


(OB-OtpJdu 

2yu D -u 


(09) 


evaluated at z = u. 

The term (<j^-0fj>)/2 in equations (C9) and (CIO) is to be 
interpreted as including • if the integration from 0 to u 3 (v) 
includes a portion of the area Sj^g. 

In the immediate vicinity of the edge u^v), the funotion A-^ 

has a value due mainly to the strength of. the pole at that- edge. 

As meas ured a long a line v = constant, in the vicinity of the 
edge ,/u-Up in equation (C9) pan be replaced by its mean value, 

in which case equation (C9) reduces to the form 
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« L. r . .. 

jr/uHigtv) 


pUgCv) : 

((Jg-Oij,)du 

' 2 v /u 3 (y)-u 


(CIO) 


&>. all cases, except for an airfoil with a pointed tip having both 
edges swept “behind the Mach' angle, the value of the integral 


^(v) 


(Cg-OgOdu ‘ 


Si (u D ) 


(Oil) 


(for example, see equation (9)), or the value of the equivalent 
expression in v will be known from, a previous calculation of 
Si(u D ) or g 2 (v). In the vicinity of the curve u *= u 3 (v), 

equation (CIO) may therefore he written as 


8i<» 3 ) 


A ~ 

1 «yu^ 


From equations (C8) and (C12) 


(C12)' 


A|S U T ? 


A « llBl fV V _ IT 

1 v-»Vg(u) * V u ~ u 3 ~ n 
Ihe contribution to gg D (v D ) from A 1 is 

r s M 


A-j_ dv 


/(v 3 -v)(v D -v) 


- A 1 log e 


V^D 7*3 


hence, gg ^(v-jj) may he written as 


S 2,7 V 7 = 


(CIS) 


(C14) 


k - -&=) -A x log, 


+ sr 


6 ^-v^3 


(CIS) 
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Finally, gg(v D ) is given "by the relations 

s 2 ( v d) -tarVa 

a) 

' ' [(g B -a;)-(a E - g a?.)2]' 


rvg(^) 


- Aj. log, 


y^ -v^ 




'3. 


dv 


- Tv „ 

: • ^ 


(0X6) 


where 


Vi 


■*/ Su V 3 


Al * Sx^) 


(C13) 


is evaluated at curve 3, 


For a purely numerical integration of the equations for Tv, 
it is desirable to eliminate the differentiation outside the 
integral sign in equation (Cl). The diaphraga slope X 2 can he 

computed by the use of the equation 


pv 


= i |r 


dv. 


D 




g(V 


rv 


1 8 (v 2 ) 





( 017 ) 



■which may be derived by integration, by parts ."prior tQ the ' 
differentiation. The 'quantity* • gg 1 * (vg)T iff flef feed, as : 


82* ( v p) = >” 82^) 


Equation (C17) is valid if g’-Cvp) is finite within the 

interval v g (u) to v and at its limit's 'and has at most a finite 
number of discontinuities..' In order to use equation (C17), the 
portions of g* (v) that do not satisfy these -conditions must be 
isolated and evaluated separately. 

The integrand of equation (C17) is always singular at the. 
point v « v D . , ' •’■■■■' " 

In general ’82* (f-n)' ’also has a pole of order l/Z at 
v = v 2 (u) contributed by the term -Kos-O^ V v D“ v 2 *4 equa- 
tion (C16) . Its' contribution to' X2 i 8 

* r ' ! 

-|(a B -a T ) 2 (ci8) 

The integral in the general expression for gg* (v D ')‘ due to 
the integral 


v^3 


- ( p b -<%)2 

2/vp^V 


■dv 


in equation (016) will have an infinite integrand at 
unless the term 



(C19) 


v D = v 2 (u). 


(C20) 


The expression (C19) may be written 
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= (^“T2) - ' av Z \ 0 


>v 2 ■; 


+ . , £aL £££*? ) dv - 


1 -'(021) 


The first terto’af ‘expression - . (C21) hue an everywhere finite deriva- 
tive with a finite integrand given by 


V '2 r^ar.-On)! 

■ (o B ^ r )-(o- B - % ) 2 - (v-T 2 ) [ - | 


4(v d -v)‘ 


dv ’ (C22) 


The contribution of the- second- term of expression' (021) .to 
s(v D ) . is 

s lr^~L': - /v^) 


• - <V v 3) A^3'] ' 


(023) 


for which the corresponding contribution'-, to gg’^p) is 




.] ~ (7 D- T 5^ + ( v D-~ 

*2 L ■ 


- J^lFz 


Plnaliy, * 2 may he computed by adding the contributions of the 
regular and singular terms of gg(vjj) and gg'Cvp): 
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>.z n i 
:v. 


rv ■ 




..... 7 .. JY 2 (h) 


•r- . .'X’ O^-Wv) 


. + ~^(v) 

where ■ -• 

. ' V' M:,’ 

1 ,r ’i ^B*^ 2 - V Td • 

h' (v d )=7“xi= + At r 




l v D -y 3 


(C24) 


r . . *J 




tv^r> : v Cvxj-v 2 );- • 




-Air^ 


K<%- 


dv J 


^ V 3 


X, - 


dy 


1 7W 


( v D -v)‘ 




2V ~- • P(ob-<j t ) 

■ (y^Hyg,)., - < T - T 2>L~'37'd 3 

•, 4(t b -t) 3 ^ . 


dy 


y,(u) 


3 . (C25) 

The function h r (vp) contains all. the regular terms of g 2 * (vp) . 

Calculation of surface yelooitles,. - The velocity potential cp 
is. glyen by the double integral . 


cp = - 5- 

^ it 


■Q d| dn- 


y(s-£) 2 -B 2 (r-n) 2 " 


(C26) 



32 


MCA IN No. 1484 


in which the integration ia .carried out over the forward Mach oone 
or over . limited .regions of ' this cone as described in reference (l) 
and where . a ’ represents the slope of the wing surf ace or the elope 
of the diaphraga, .as required.’ ’ : 


In most cases interest centers in the velocities or in the 
pressure coefficients, which are proportional to the x component 
of the velooity. She components of velocity are given by the 
equations 


V 

* " Sr ~ 


ui 

* dx 


y(i-i) 2 -P 2 (y-ri) ! 


Sep 

^ * Sr 


E iL 

*. Sy 


od| dn 


//(2-.|)^8 2 .(y-Ti) 2 


(C27) 


The manipulation of these derivatives is facilitated by transform- 
ing the integrals to the u,v coordinates . •' 


>u w 




- - 


tr 

Mjf 


JO J< 


gdu dv 

}(v w -v) 


(C28) 


inhere u^’ and v w are the . coordinates of a point on the wing. In 
terms of equations (C27). arid (028) 


*-§?-&(!£♦!£)■ 

Y = £\ 
y dy 2 \ov w ou^j. J 


(C29) 


The derivatives of the potential with respect to 
are given by the relations’ • 


% and % 



NACA TN No. 1484. 


33 ■ 


kU w 


rv, 


w 


Sep -g S 
= Su^ 


a du dv 


*/ ( v u) /v*' 


nu^ rv,. 


Sep ^ '-g _S_ 
Sy w C Mjt ^r w 


JO J( 


a du dv 


*/%-u V (%-v) 


(030) 


The calculation of the derivatives may he simplified in -the 
maimer subsequently decrihad for the calculation of Scp/S v^ T , The 
equation for the potential may he Witten as 


pttw 


Pv 


qp 



du 

</v^ 


^idv 


(C3i) 


where the symbol v w is tmderstood to include v-p when the limit 

of integration lies on one of the diaphragms, in whioh case the 
contribution of o from both top and bottom wing surfaces forward 
of Vp must he considered. 


It is convenient to define a function B(u,v^) by the 
relation 


B(u, 7 w ) 


’w 


q(u,v) dv 


r V -V 

¥ W v 


(C32) 


then 


St? v_ 

Su^. 2 Mjt 


>u w 


B(u,v w ) - E(%,v w ) E(%,v v ) 

— du - 2 






(C33) 


For points on the top giurfaoe of .the wing p^u^Yy.) is given 
by the relation 
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*= 2<3j(u,V tf ) y V w rV^ + A^Cu) logg ~J| 

*/ v w “ v v 3 



dv 


,vV y 


‘ <%(^v) - cr T (u,v w ) ^ 

+ — tt 7-- — dv 


and the value of Rg(u,v w ) is gijea by the relation 


e b( u * t w ) *■ 2 cj b(^ v w ) - A 1 (u) log Q 




v v " , /Tj 


*1 - 


A l (u) ] dv 


/v 3" v y v w - v 


w 

+ - ofe(u,v v ) 


For a point on the after diaphragm (A g ), the value of 
B g?( u >' v I)) s*“L ”%(^Td) 18 

= - %(u,v D ) « -2(03+03, ^/v D -v 5 - yvp-vg ^ 
■'Vg(u) 

* + - [fo B +qp) - (a F KT T ) 2 ] ~=z 


fr V-v 
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For the calculation of lift distributions, the contrihutibn of 
%( u /Td) 01131 Eg(u,Vu) cancel ‘and- may therefore he neglected. 

For flat -plate wings. Cg + arp » 0. 

For points on the forward diaphragn (A-^), the .-defining 

integrals for ECu^Vp) includes only the effect of (Che 

value of : 'E 'is therefore zero'. • 

The Integrand of the integral in equation (C33) generally 
becomes infinite at u = -u v * Although -this integral- could be . y 
integrated -by parts and the’ integral obtained in terms of the" 
derivative of E(u w; y w ) with respect to 'it Beeins preferable 

in this 'case to expand E(u,v v ) in a power series about the : 
point u^., "When the function B(u-,v^) can be expanded with suf- 
ficient accuracy in a segment of width 25 _ 

>‘u.5i u-25 - ‘ ’ ’ ■ 

* ' « » 
as a parabola of the form . . - . . . >. 

E(u,v v ) = E(%,v v ) + Aiu^-u) + B(%-*) 8 

and if 

P(u, r -25,v w ) - B(iV,v w ) = a 
R(%-8,v w ) - B(u v ,v w ) = b 

it may be easily shown that 

R(^\.)-R(iV,v w ) / 2 " gb -a 

K-u)3/B 3 



It is of interest to note that the function B is finite for 
all values of u except u = U 3 (v). B(ug,v w ) and are 

OXLy 

finite at the point (u 3 ,v w ) if A 1 (u 3 ,v 3 ) is zero. 

The calculation of ^2. is carried out similarly to that of 
dcp ^ T w 

^-2-, except that the roles of u and v are interchanged and a 
suitable first integral is defined to take the place of B(u,v w ). 
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Fig. 6 
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Figure 6. - Wing plan form and diaphragm areas for equations (17) and (18). 
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Figure 0. - Concluded. Diaphragm slopes Va for symmetrical flat-plate delta wing. 
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